PROCEEDINGS 

OF SCIENCE 



o 

O 

Q 

(N 



> 

(N 



X 



Exploring Gauge-Invariant Vacuum Wave 
Functionals for Yang-Mills Theory 



Hilmar Forkel* 

Institut fiir Physik, Humboldt-Universitat zu Berlin, D- 12489 Berlin, Germany 
E-mail: 



f orkel@physik . hu-berlin . de 



We study gauge-invariant approximations to the Yang-Mills vacuum wave functional in which 
asymptotic freedom and a detailed description of the infrared dynamics are encoded through 
squeezed core states. After variationally optimizing these trial functionals, dimensional transmu- 
tation, gluon condensation and a dynamical mass gap of the expected magnitude emerge trans- 
parently. The dispersion properties of the soft gauge modes are modified by higher-gradient 
interactions and suggest a negative differential color resistance of the Yang-Mills vacuum. Cast- 
ing the soft-mode dynamics into the form of an effective action for gauge-invariant collective 
fields, furthermore, allows to identify novel infrared degrees of freedom. The latter are gauge- 
invariant saddle-point fields which summarize dominant and universal contributions from various 
gauge-field orbits to all amplitudes. Their analysis provides new insights into how the vacuum 
gluon fields generate gauge-invariant excitations. Examples include a dynamical size stabilization 
mechanism for instantons and merons, a gauge-invariant representation of their effects as well as 
a new physical interpretation for Faddeev-Niemi knots. 



International Workshop on QCD Green 's Functions, Confinement and Phenomenology 
5-9 September 2011 
Trent o, Italy 



* Speaker. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http://pos.sissa.it/ 



Exploring Gauge-Invariant Vacuum Wave Functionals for Yang-Mills Theory 



Hilmar Forkel 



1. Introduction 

The Hamiltonian formulation of Yang-Mills (YM) theory in the Schrodinger picture, although 
not particularly efficient in the perturbative domain, offers considerable benefits when addressing 
nonperturbative issues. Among its attractive features are the explicit representation of the vacuum 
state which invokes quantum mechanical intuition [|l|], the ability to treat genuine real-time prob- 
lems (including non-equilibrium processes) as well as the transparent description of topological 
effects [f2|]. In particular, however, it makes gauge theories accessible to a variational treatment 
[jl|, [3p, i.e. to one of the few approximation schemes currently available for strongly coupled quan- 
tum field theories. 

Variational calculations in Yang-Mills theories are often performed in a fixed gauge, most no- 
tably in Coulomb gauge [Q]. In the following we will report on our complementary explorations 
[]|] of a manifestly gauge-invariant formulation of the variational problem [§J. This framework 
renders fundamental infrared (IR) physics, including dimensional transmutation and the genera- 
tion of a mass gap, particularly transparent. Moreover, it preserves the full topological structure 
of the gauge group. The latter is particularly relevant since topological properties are likely ro- 
bust enough to survive limitations of the restricted trial functional basis which keeps the approach 
analytically manageable. Another attractive feature of the gauge-invariant formulation is that the 
infrared dynamics can be re-expressed in terms of gauge-invariant collective fields which subsume 
contributions from whole gauge-field orbit families. After performing an IR improved variational 
analysis we will make use of this feature to identify gauge-invariant and universal IR degrees 
of freedom of the gauge dynamics [Q]. More details can be found in Refs. [§, 0], 



2. Gauge-invariant vacuum wave functionals 



Starting from an approximate and hence typically gauge-dependent "core" functional Yo 
of the static gauge fields (i.e. of half of the canonical variables), we impose gauge invariance by 
averaging over the gauge group. The result is a trial vacuum wave functional (VWF) of the form 

(2.1) 

where d/j. is the Haar measure of the S\J(N C ) gauge group, Q the topological (homotopy) charge of 
the group element U&\ and 6 the vacuum angle. Since the vacuum wave functional is nodeless [fit], 
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and expand the real functional <I> into a power series. 
The constant term is absorbed into jV and the term linear in A is generally discarded (coherent 
gluon vacuum states are known to be unstable [§]). 

The next term is quadratic in A and plays several crucial roles. First, it removes the ambiguity 
in *Po [||] due to the invariance of the Haar measure in Eq. (2T) under group transformations. 



Furthermore, this term can incorporate asymptotic freedom and thus render the VWF exact in the 
ultraviolet. Finally, and from the practical perspective most importantly, functionals resulting from 
a quadratic term can be integrated over A analytically. Hence one generally truncates the series for 
<1> after the quadratic term, which leads to the "squeezed" core functional 
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with the normalization factor jV g 1 



[det (G/2)] ^ 4 and a real "covariance" G 



-l 



2.1 Gluon dispersion: asymptotic freedom and IR generality 

We now have to specify the properties of the function G 1 in the trial functional family (2.2). 



Translational invariance was already anticipated in Eq. (|2.2|). We will further restrict ourselves to 
a purely transverse covariance with the Fourier transform & 



G, 



-l.ab 
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(2.3) 



(cf. Ref. [Q] for a discussion of this choice and Refs. [ JTOj , [TT| ] for the impact of longitudinal 
contributions). The normalizability of physical wave functionals then demands G 1 (k) > and 
further ensures vacuum stability and a positive energy spectrum. In order to implement the correct 

' by 



(G. 



Vo 



(G 



UV behavior, we factorize the core functionals ( |2.2| ) as 

splitting the k integration domain in their exponentials into soft/hard momentum regions with 
jU. The separation scale \i will be determined below. Asymptotic freedom requires G to approach 
the non-interacting, massless static vector propagator Go(k) = l/k for k — > oo. As long as jJ. S> 
Aym (where Aym is the Yang-Mills scale) perturbative hard-mode corrections remain small, which 
allows us to approximate 

G- 1 (k) = G Q 1 (Jfc) = k. (2.4) 

The unknown IR covariance G<* (k), on the other hand, will be determined variationally. In Ref. 
[^] the minimal one-parameter trial function G~ l KK (k) = \i was adopted. We have implemented a 
far more comprehensive parametrization based on the under reasonable analyticity assumptions 
general and controlled gradient expansion 



G-\x-y) 
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(2.5) 



Eq. (2.5) can be efficiently truncated to maintain an analytically manageable trial basis for the 



soft-mode physics. Besides jj,, the variational parameter space now contains the IR gluon mass 
m g > and a few low-momentum constants c; which characterize dispersive gluon properties in 
the vacuum. The regularized delta function 5< (x—y) := / d 3 k/ (2k) 3 6 (jl 2 — krj e^~^ encodes 
the slow variation ||dA<|| / ||A<|| < jJ. of the soft modes and ensures that the higher-order terms in 



Eq. ( J2.5[ ) are parametrically suppressed. 

As a consequence of G 1 (k) > 0, the low-momentum constants are subject to the bounds 
c\ < 1, c% > — 1, etc. (for m g > 0). Requiring continuity of G 1 (k) at the matching point k = }X, 
furthermore, fixes m g as a function of the other variational parameters. When truncating to c,->2 = 0, 
for example, one has 

m g (p,ci) = - . (2.6) 

i — c\ 

Note that the requirement of a non-negative IR gluon mass restricts the c\ domain to c\ < 1, in 
agreement with the above bound from G 1 (it) > 0. The VWF ( |Q| ) together with the core func- 
tional (2.2) and the covariance (2.4), (|2~3|) (possibly with perturbative corrections) appears to be 



the "richest" gauge-invariant trial functional family whose matrix elements can be calculated ana- 
lytically by currently available techniques. 
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3. Variational analysis 



One the basis of the trial functional family (2.1) discussed above, the variational analysis 
amounts to minimizing the expectation value 



(JT(A,E)) = 
of the Yang-Mills Hamiltonian density 
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(3.1) 



(3.2) 



(in temporal gauge, with E a = iS / 8A a ) with respect to the parameters appearing in G l . After 
inserting the wave functional ( |Q| ) into Eq. (3T) and interchanging the order of integration over 
fields and group elements, the gauge invariance of the A integral allows to factor out a gauge group 
volume. Eq. ( fO| ) can thus be rewritten as 

\DAy Q \A u ]j#> U a ,A) XjfoU 
{J? (A,E)) = 1 J dA r ^ L -^- (3.3) 
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(where DU is the functional S\J(N C ) measure as defined in Eq. (gjj)). After evaluating the func 
tional derivatives contained in Ji?, the Gaussian integration over A can be performed exactly, re 
suiting in 

SDU(((JiP)))ex V {-r b [U]} 



fDUexp{-r h [U}} 



(J?) 

where we introduced the notation 

[A...A...E...E } ))exp{-F h [U}} = I DAy 
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(3.4) 



(3.5) 



for matrix elements between U -rotated and unrotated core VWFs. The above expression defines, in 



Wo 



which describes dynamical 



particular, the effective bare action [U] = — In J DA 
correlations generated by the gauge projection. This action gathers all those gauge-field contribu- 
tions to the generating functional whose approximate vacua Yo at t = ±°° differ by a relative gauge 
orientation U. Hence the gauge-invariant "variable" U represents the contributions from all such 
gluon field orbits to the vacuum overlap. Explicitly, one finds [|]] 



r b [U] = ^- I Jd 3 xJd 3 yLl{ 
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G <5 +... where G t 



■ G ab {x-y) 



with Li = WdiU =: LfX a / (2i) and 2> ab = (G + G u ) 

(x) (A72) U(Z)®U (y) (X b /2) (y). 
After splitting U (x) = £/< (x) £/> (x) with U> (x) = exp (-ig(j) a (x) X a /2) into hard- and soft- 
mode contributions and integrating over the hard modes <p" perturbatively jjj, furthermore, one 
arrives at 

SDU < !D<S>(((3t?)))z W {-T b [<$>,U < ]} 



(3.7) 
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With the additional definition 



exp{-r< [t/<]}:= / 



which contains the effective soft-mode action 



exp {-r fo [0, £/<]}, 



T< [!/<]:=- In y"D0exp{-r fc [0,I/<]} 



(3.8) 



(3.9) 



(i.e. the RG evolved bare action fl3.6| )), we can finally rewrite the matrix element ( J3 - 1[ ) solely in 
terms of the t/< field dynamics, i.e. 



jD^<((jr))exp{-r< [£/<]} 
^ ; /D^<exp{-r< [t/<]} ' 



(3.10) 



Since the "reduced" (i.e., fixed f/<) matrix element {{J^}} is a nonlocal functional of the soft 



modes U < , the evaluation of Eq. ( |3.10[ ) amounts to calculating (equal-time) soft-mode correlation 
functions 



3.1 Vacuum phases 

In integrals over the U < fields (such as those in Eq. ( 3.10 )) the unitarity constraint t/<t/< = 1 
can be resolved by inserting a delta functional which is then written as an additional integral over 
a hermitean auxilary field E. In Eq. ( |3. lOj ) the integration over the then unconstrained U < becomes 
Gaussian and can be done analytically. In the mean-field approximation, the expression for ((<%?)) 

— _ 9 

is then evaluated at the saddle point E =: (jU,£) of the E integral, i.e. at the minimal-action solution 
of the gap equation 

'ul AB (x)U < , BC (x)) = 8 AC (3.11) 



which reintroduces unitarity in the mean. After adopting the one-loop Yang-Mills coupling y (ju) = 
g\ M {}x)N c /lz 1 = 24/ (IHqju/Aym)) tne solutions f of Eq. ( |3.11 ) depend on two variational 
parameters, the RG scale pL > and ci < 1. The critical line jU c (ci), i.e. the parameter subspace 
where the (dis-)order parameter f (fi c (ci) ,c\) vanishes and the phase transition takes place, can be 
found analytically as the combination of the two curves 
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(l-c 1 )r(ci)±^5r 2 (ci)(l-ci) 2 -4(l-c 1 )[l-c I F(c 1 )] 



(3.12) 



(f (ci) := arctanhy^/y^). We plot this closed phase boundary in Fig. [Ij It limits the parameter 
ranges to 0.5 < < 8.86 and —0.48 < c\ < 1 and thus prevents the minimal-energy solution 
I* from attaining unacceptably large values of pL and |ci|. Nonzero solutions of the gap equation 
exist only when the gauge coupling exceeds a critical value, i.e. for g 2 (ju) > g 2 (c\) , as expected 
on physical grounds. The (dis-)order parameter goes to zero continuously, furthermore, i.e. the 
disorder-order transition is of second order (which may be an artefact of the mean-field approxi- 
mation 
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-0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 



Figure 1: The vacuum phase diagram. Inside the plotted phase boundary jj. c (c\) /Aym the theory is in its 
strongly-coupled disordered phase. (The underlying approximations are reliable for fi > 4Aym and ci 6 
{-0.5,0.5}). 



3.2 Vacuum energy density 



Working with the Poincare-invariant trial states (J2JJ) and taking only one-loop corrections 
from the hard modes into acccount, it is sufficient to regularize Eq. ( |3.1[ ) by a momentum cutoff 
Auv [@L Separating the complete vacuum energy density £ = E/V = (Jym) into hard and soft 
contributions, 

e(ji,ci,£;|) = (JeW) = e>(M) + e<(M.ci,C;l), (3-13) 
(£ = trig/ /i) the cutoff dependence resides solely in 
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(3-14) 



As expected, this is the (regularized) zero-point energy density of two transverse, massless vector 
modes in the adjoint representation of SU(/V f ) with energy a> (k) = k. Simple normal-ordering thus 
subtracts the Auv dependent term. For c,->2 = one then finds the total energy density e (jU,ci) := 
£ (H,ci, £ ct (ci) (ju,ci)) in the disordered phase as 

"4c^ + 10c? - 50ci +30 1 - ci arctanhy^cT 



A? 2 

e(M,ci)=-^M 



30c 2 (l-ci) 



+ 



? 2 - 2ci? 3 + c 2 f 4 + 2yc t (1 - ci) ? 2 (j3 - 2ci74 + c 2 j 5 ) 

i 3 ^ 



(3.15) 



where the integrals f n (£,ci), j n (£,ci) are defined in Ref. [^] and evaluated at f (jU,ci). This 
energy density is plotted in Fig. ||. 

In the ordered phase, i.e. for /i ;g> Aym where g 2 (/j)<Cl, the energy density can be calculated 
perturbatively (in g ). Since fluctuations <p" around U < ~ 1 are small in this phase, one may 
approximate U < = exp (igcp^X") = 1 + igcp^X" + (g 2 ) . After adding the hard-mode contribution 



( |3.14[ ) and discarding the zero-point contribution, this results in 



£(jU,Ci 



Nr. 



An 2 



1 4 1 



c\ + 15c 2 — 50ci +30 arctanh-y/cT 
30(1 -ci) 2 



(3.16) 
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4 




-0.5 0.0 0.5 



Figure 2: The energy density e (jlt,ci) of the vacuum field solution \ (n,c\) in the disordered phase. (The 
plot shows the parameter ranges jx g {4,9} Aym and c\ e {—0.5,0.8} .) Note the minimum of the energy 
surface at c\ ~ 0.15. 



It is reasonable to expect that this perturbative result remains qualitatively reliable down to the 
phase transition at jU c [0]. The singularity of the energy density at c\ — > 1 encodes the vacuum 
instability for c\ > 1 and thus automatically ensures that the wave functional remains normalizable 
during the variational analysis. 

The most important lesson of the above analysis is that £ (jU,ci) increases monotonically with 



/i and c\ (for — 2 < ci < l)in the ordered phase while the energy density (3.15) in the strongly- 
coupled disordered phase monotonically decreases with /I and ci, up to the phase transition. This 
indicates that the vacuum energy density becomes minimal at the phase boundary in the disordered 
phase, i.e. at f = + (where the number of massless particles becomes maximal ^). The pre- 
cise minimum, e{}X* ,c\) ~ — 210.59A YM , is reached at c\ ~ 0.15 with /I* = ll c {c\) = 8.61A Y m- 
These values justify the perturbative treatment of the hard modes and of the 4U contributions. The 
c\ corrections reduce the vacuum energy density by about 11% and provide a rather substantial 
improvement of the wave functional. 

3.3 Gluon condensate and quasigluon kinetic mass 

At the physical parameter values, i.e. at the border of the disordered phase where the energy 
is minimal, the gluon condensate becomes 



, 2V Af?-1 4 



lc\ - 20fc\ + \5c\ + 20fc\ - 50ci +30 ( 1 - c x 2f\ arctanh^/Q 



30c?(l-ci) \ c\ 3 J y/Ti 

(3.17) 



(7* = g 1 (n*)N c /n z ~ 1.012). Numerically, this implies 

—F 2 \ = 20.87Ay M ~ 0.01 1 GeV 4 (3.18) 

(for Aym — 0. 15 GeV), i.e. an about 25% larger value than in the uncorrected c\ = case. The re- 
sult ( p.!8| ) lies comfortably within the standard range {(a/n)F 2 ) = 0.0080-0.024 GeV 4 obtained 



from QCD sum rules [12] 
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Our finite and positive result for c\ has further interesting consequences since it reshapes the 
composition and dispersion of the vacuum field population. Indeed, the attractive IR interactions 
generated by c\ > deplete the density of ultralong-wavelength k — > modes and populate the 
k ~ ji modes more strongly. This is consistent with the expected average wavelength X ~ Ay^ of 
the vacuum fields. Since G 1 (k) describes the dispersion relation ft) (k) of "quasigluon" modes in 
the vacuum, furthermore, one may relate c\ to the modulus of the dimensionless quasigluon group 
velocity v(t\ = dG^ 1 (tj /dk at k = jj, [§], 

I I v (a*) niQ , 
ci = tttn i v ( 3 - 19 ) 
v(ju) + 2 

For > c\ > 1 (as in our case), furthermore, the "effective kinetic gluon mass" m g , which relates 
velocity and momentum as k = m g v, is negative. Hence v is opposite to the momentum, causing 
the "quasigluons" in the vacuum to decelerate when an external force is applied. (Such dispersions 
are encountered in several condensed-matter systems and are in stark contrast to the behavior of 
free gluons.) Hence quasigluons (with their small scattering amplitudes) may show a negative 
differential color resistance. 



4. Infrared degrees of freedom 

Our above representation of the vacuum dynamics in terms of the gauge-invariant low-energy 
fields U < provides the opportunity to search for specific U < which may play a particularly impor- 
tant or even dominant role in the generating functional (and hence universally in all low-energy 



amplitudes) [13]. If such fields exist, they can be regarded as universal infrared degrees of freedom 



(IRdofs). In contrast to other proposed IRdof candidates (e.g. classical gauge-field solutions like 



instantons [|14[], or monopole and vortex configurations), the IRdofs expressed in terms of U < are 
gauge invariant and contain crucial quantum effects (e.g. those which stabilize the instanton size, 
see below). From a practical perspective, these IRdofs will be useful as well since many technical 
problems encountered when dealing with gauge-dependent fields are avoided from the outset. Be- 
low we will show that large classes of such IRdofs indeed exist and review how their stability and 
topology emerges. We then construct important IRdof classes explicitly and discuss their properties 
and physical interpretation. 

4.1 Gauge-invariant saddle point expansion 

We start from the vacuum overlap matrix element, i.e. the functional integral 

Z = yDt/ < exp(-r[t/ < ]) (4.1) 



over the soft modes, with the action T given by Eq. ( ]3.9| ). (Sources can be included when needed.) 
A steepest descent approximation for Z can be set up by expanding U < around the saddle point 



fields Uj (x), i.e. the local minima of the soft-mode action (3.9) which solve 



5r[£/<] 
5<y< (x) 



u<=u. 



IQ) 



■ 0. (4.2) 
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(Different topological charges (see below) are summarily denoted by Q since the action is varied 
in each topological sector separately.) To leading order, the saddle point expansion for Z is then a 
weighted sum (or integral - the symbolic label i becomes continuous when the saddle points form 
continuous families) over the contributions from all relevant solutions 



L F i 

QeZj 



exp 



(4.3) 



where nontrivial pre-exponential factors f) are typically generated by zero-mode contributions. 



For the general analysis and explicit solution of Eq. (i.2) we adopt the parametrization 
U < (x) = exp [0 (x) h a (x)X a / (2i)] of the S\](N C ) elements and work directly with the N 2 - 1 in- 
dependent degrees of freedom of U < , i.e. the unit vector field h a and the spin-0 field 0. For 
simplicity, we will also specialize to N c = 2 and use the first two terms in the expansion ( |2.5| ) of 
the inverse finite-mass gluon propagator G 1 (k) = yjk 2 + jx 2 as a template for the covariance [0]. 
The soft-mode Lagrangian can can then be written as a sum of two- and four-derivative terms, 

ja?(tf<) = jSf M (M)+^i(M). (4-4) 

(For the explicit expressions see Ref. The saddle point equation (fO|), when specialized to 
variations with respect to <p and n" , becomes a system of four nonlinear partial differential equa- 
tions. Its localized solutions can be shown to be stable under scale transformations, due to the virial 
theorem (1) = (1) [01. (Clearly the four-derivative term is crucial here - truncation of 
the gradient expansion ( |2.5| ) to two powers of dU < /jl is therefore the minimal approximation 
which supports stable saddle points.) The origin of this stability can be traced to the mass scale pt 
emerging from the out-integrated short-wavelength quantum fluctuations. 



An already mentioned, crucial benefit of the gauge-projected wave functionals (|2JJ) is that 
they fully implement the nontrivial topology of the gauge group and fields. The U < fields thereby 
inherit three integer topological quantum numbers [f7|]: a winding number 6[t/<] (characterizing 
the homotopy group 713 (S 3 ) = Z), a monopole-type degree q m [n] based on 712 {S 2 ) = Z and finally a 
linking number qn [n] in the Hopf bundle Kj, (S 2 ) = Z which classifies knot solutions. This topology 
entails two lower action bounds []7|] of Bogomol'nyi type, 

r[t/<] > * leifcji, r[fc = (2*+i)M] > 2r „, \i*m v \ (4-5) 
g W g w 

which ensure that contributions to soft amplitudes from saddle points in high charge sectors can 
generally be neglected. This allows for practicable truncations of the saddle-point expansion. 
(Saturation of the first bound requires the fields to solve the Bogomol'nyi-type equation djLj = 
^pLEijkLk, incidentally, which can be considered as the analog of the self-(anti)-duality equation in 
Yang-Mills theory.) 

4.2 Important examples of gauge-invariant infrared degrees of freedom 

In general, the saddle-point solutions have to be found numerically. Among the exceptions are 
the translationally invariant vacuum solutions U c = const, (which are the absolute action minima 
r [U c ] = 0) and several nontrivial solution classes which can be found analytically. In addition, 
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important and sufficiently symmetric solutions classes can often be obtained by solving substan- 
tially simplified field equations [0]. (The typically smaller action values of solutions with higher 
symmetry generate a stronger impact on the matrix elements, furthermore.) 

As an example for nontrivial analytical solutions, we consider U < fields with constant n a for 
which the saddle point equation becomes linear: d 2 (d 2 — 2/i 2 0) = 0. The general solution does 
not carry any topological charge and was found in Ref. The subset of spherically symmetric 
solutions with finite action, in particular, is 

^n=c) {r)=Ci + ^^( l _ e -^r\ (46) 

\J2\Xr V / 

with the action T [^ h=c \n c ] = c\%j (V2g 2 (/in. Since Eq. ( |4.6| ) is not subject to topological 
bounds, it continuously turns into one of the vacuum solutions for C2 — > 0. 

A particularly important saddle-point solution class consists of topological solitons of "hedge- 
hog" type, 

n a (x)=x a , (f, (x) = (f> (hh) (r) (4.7) 

(x" =x/r, r = \x\). Well-defined hedgehog fields must satisfy the boundary condition <j)^ hh > (0) = 
2ki7i (regularity at the origin further requires 0" (0) = 0) and finite-action fields additionally have 
(f>( hh ) (oo) = 2&2 ft where and the charge Q = k\ — ki are integers. The more general boundary 
conditions 



dhh) 



n — m 



(4.8) 



(n,m integer) additionally admit infinite-action solutions with half-integer winding numbers Q (for 
either m or n odd). All hedgehog fields further carry the monopole-type charge :=q m [x] = ±l. 
Due to the periodicity in (j), it is sufficient to consider boundary values in the range (j) (0) G ]0,2tt]. 
The dynamics of (f) (r) is governed by the radial Lagrangian 



g 2 0)M 



X - [rf) 2 + (3 + At 2 r 2 ) (f ) 2 + 4.U 2 ( 1 - cos 0) 



(4.9) 



The hedgehog saddle points, found numerically in Ref. [0], turn out to comprise mainly contribu- 
tions from and around the gauge orbits of the classical Yang-Mills solutions, i.e. (multi-) instantons 
and (multi-) merons. The potential term in Eq. (|^) is analogous to that of a one-dimensional pen- 
dulum in a gravitational field, with stable (unstable) equilibrium positions at = % ((j) = 0), modulo 
multiples of 271. 

We first discuss the regular hedgehog solutions. Their three boundary conditions (0) = 27T, 
<p" (0) = and <p (oo) = 2% (1 — Q) imply that for a given Q all of them can be found by varying 
the initial slope j8 := 0'(O). (For the irregular solutions with (0) = K see Ref. [0].) The reg- 
ular solutions turn out to contain one. finite-action solution for each Q, denoted as the "\Q\ (anti) 
instanton class", and the remaining, continuous (in j3) infinite-action families, the "2\Q\ (anti) 
meron classes". The 1 -instanton class solution is depicted in Fig. |3[ Its relative gauge orienta- 
tion U = UZ 1 U + is even quantitatively close to that of the Yang -Mills instanton [14] (orbit) (see 
dashed curve in Fig. ||). This confirms that the \Q\ instanton classes indeed primarily summa- 
rize Yang-Mills instanton contributions. However, they also contain crucial, dilatation-breaking 
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Instanlon-class solution (0=1) 


















♦ liVM (r)=2n[1-r/(t'-4n 
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Figure 3: The 1-instanton class solution. The dashed line corresponds to the Yang-Mills instanton. 



quantum corrections which dynamically stabilize the instanton size at about p ~ 2\i 1 , compat- 



ible with instanton liquid model [ |14[ ] and lattice Q15Q results, and thus overcome the chronic IR 
instabilities of classical Yang-Mills instanton gases. Since instanton effects play important roles 



in Yang-Mills theory (e.g. in the 6 vacuum [16] and in spin-0 glueball physics [12, tt4, 17]), it is 



crucial that they are (at least partly) included in the vacuum functional ( |2.1[ ). In fact, approximate 
non-hedgehog solutions corresponding dominantly to ensembles of instantons and anti-instantons 
should also exist and play prominent roles (since they would be enhanced by a large "entropy", as 



in phenomenologically successful "instanton liquid" models [|14]]). 

All remaining regular (i.e. <p (0) = 2jt) hedgehog solutions, with initial slopes j3 between the 
discrete instanton-class values p\g, form the 2\Q\ (anti-) meron classes. Those approach one of 
the values <pM (°°) = (2& + 1) 71 at spacial infinity and therefore have infinite action, as the Yang- 
Mills merons [|T6|]. Moreover, solutions with <p (°°) = (2k + 1) jt, corresponding to an odd number 
of merons, carry the half-integer topological charge Q of their Yang-Mills meron counterparts. In 
addition, quantum effects ensure that our meron-class solutions acquire a finite size and therefore 
remain nonsingular. Since the 2|<2|-meron classes appear in continuous families (parametrized 
by their "size" j8 _1 ), furthermore, their large entropy will help to overcome their infinite-action 
suppression in functional integrals. As in Yang-Mills theory, merons could then play a physical 



role, e.g. in the confinement mechanism [16]. 

We conclude our discussion of selected saddle-point solutions with one of the most intriguing 
classes, consisting of (solitonic) links and knots. Those emerge from a generalization of Faddeev- 
Niemi theory Q], 

{diff f + \ (e ahc dih b djff^ 2 + [s abc n h d 2 n c ^ 2 



; (a0 



(4.10) 



which is included in our soft-mode Lagrangian for constant 0& = (2k + 1) jr. The corresponding 
solution classes U < (n) describe twists, linked loops and knots made of closed color fluxtubes. 



Since Eq. (1.10) follows uniquely from the VWF (2.1) and the Yang-Mills dynamics, our approach 
provides a new framework and physical interpretation for such solutions. In fact, they remerge 
as gauge-invariant IR degrees of freedom representing sets of gauge-field orbits with a collective 
Hopf charge. While the n field of the Faddeev-Niemi model is interpreted as a gauge-dependent 
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local color direction in the vacuum, in particular, our n is manifestly gauge-invariant. This may put 
the tentative interpretation of such knot solutions as glueballs on a firmer basis. 

5. Summary and conclusions 

We have studied gauge-invariant wave functionals for the Yang-Mills vacuum which incorpo- 
rate asymptotic freedom and an a priori general dispersion for the infrared gluons in Gaussian core 
functionals. In this at present probably richest analytically manageable and gauge-invariant trial 
functional basis, we have then variationally determined several vacuum properties. Dimensional 
transmutation, dynamical mass generation and gluon condensation emerge transparently and gen- 
erate mass scales consistent with other approaches. In addition, the improved vacuum description 
in the infrared predicts a negative kinetic mass of the soft gauge-field modes and thus suggests a 
negative differential color resistance of the Yang-Mills vacuum. 

Another benefit of the gauge-invariant framework is that the dynamics can be reformulated as 
an effective theory which represents sets of gluon orbits as gauge-invariant matrix fields subject 
to higher-gradient interactions. In this effective theory we have set up a saddle-point expansion to 
determine the collective fields with maximal impact on functional integrals. These saddle points 
play the role of gauge-invariant infrared degrees of freedom. They are stabilized by the dynamical 
mass generation mechanism and inherit a rich topological structure (three topological charges and 
action bounds of Bogomol'nyi type) from the Yang-Mills gauge group. Moreover, they provide 
the principal input for a systematic saddle-point expansion of soft amplitudes (such as glueball 
correlators). 

Several of the more symmetric and important saddle-point solution classes have been found 
explicitly. Among them are topological solitons related to the classical Yang-Mills (multi-) pseu- 
doparticle solutions which mediate tunneling processes in the vacuum. Those generate a gauge- 
invariant representation of instanton and meron effects which includes quantum fluctuations. The 
latter stabilize the pseudoparticle sizes dynamically, in particular, and thereby cure the notorious 
infrared deseases encountered in dilute instanton ensembles. Similar solutions with other types of 
topological charges exist as well but seem to have no obvious counterparts in classical Yang-Mills 
theory. One of the most intriguing solution classes, finally, consists of solitonic links and knots. 
Those emerge from a generalization of Faddeev-Niemi theory which turns out to be embedded in 
our soft-mode dynamics. Hence in our framework the knot solutions find a new and in particular 
gauge-invariant physical interpretation, potentially related to glueballs. 
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